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Classical  solid  oxide  fuel  cell  anode  (Ni-cermet)  could  be  employed  as  solid  oxide  electrolysis  cell  cath¬ 
ode.  Ni-cermet  has  been  synthesized  and  tested  as  solid  oxide  electrolyzer  cathode  using  three-electrode 
techniques  between  700  °C  and  900  °C.  yttria  stabilized  zirconia  was  used  as  the  electrolyte  and  Pt  as  the 
counter  electrode.  Polarization  curves  and  impedance  spectra  have  been  analyzed  for  two  gas  com¬ 
positions.  The  presented  results  demonstrated  an  influence  of  Ni-cermet  electrode  behavior  upon  gas 
composition  and  temperature.  The  present  results  highlight  a  mechanism  changing  on  Ni-cermet  elec¬ 
trode  upon  gas  composition.  In  a  second  part,  a  one-dimensional  steady  state  model  is  developed  to 
predict  the  cathodic  behavior  of  Ni-cermet.  This  model  takes  into  account  mass  and  charge  conserva¬ 
tion,  transport  of  species  and  reaction  kinetics.  It  considers  the  porous  electrode  to  be  a  homogeneous 
medium  characterized.  The  influence  of  varying  chemical  and  electrochemical  steps  kinetic  on  the  shape 
of  polarization  curves  is  discussed.  At  high  overpotential  values  the  model  with  two  rate-limiting  steps 
has  been  validated  using  numerical  optimization  method. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Significant  research  efforts  are  achieved  to  develop  hydrogen 
economy.  Nowadays,  industrial  mass  production  of  hydrogen  is 
mainly  based  on  hydrocarbons  reforming.  However  water  electrol¬ 
ysis  could  be  the  most  convenient  production  process  if  it  uses 
a  clean  renewable  energy  source.  Solid  oxide  fuel  cells  (SOFC) 
have  been  intensively  investigated  during  this  last  decade.  Major 
improvements  have  been  achieved  and  SOFCs  now  demonstrate 
high  performance  [1],  Solid  oxide  electrolysis  cell  (SOEC)  has  been 
studied  since  1980s  as  a  promising  way  to  produce  massively 
hydrogen  [2],  High  temperature  operation  appeared  more  efficient 
due  to  the  decreased  open  circuit  voltage  and  polarization.  Using 
this  technology,  a  low  hydrogen  production  cost  is  aimed  and  a 
recent  study  indicates  that  SOEC  technology  has  a  promising  poten¬ 
tial  for  hydrogen  production  from  renewable  energy  sources  [3]. 
The  cell  consists  of  the  assembly  of  a  three-layer  region  involving 
two  ceramic  electrodes  separated  by  a  dense  ceramic  electrolyte 
made  in  the  same  materials  as  for  a  SOFC.  The  hydrogen  electrode 
is  usually  composed  of  nickel  and  yttria-stabilized-zirconia  (YSZ) 
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cermet.  The  electrolyte  is  made  of  YSZ  and  the  oxygen  electrode 
is  based  on  perovskite-type  oxides,  which  is  usually  strontium- 
doped  lanthanum  manganite  Laj.xSrxMnOs  (LSM).  The  properties 
of  these  materials  were  extensively  studied  and  are  reasonably 
known  today.  However  several  problems  are  still  unsolved:  the 
electrical  and  electrochemical  stability  of  materials  (electrolytes 
and  interconnects)  in  both  oxidizing  and  reducing  atmospheres 
and  the  mechanical  stress  due  to  thermal  behavior  of  SOEC.  More 
interest  arises  for  SOEC  with  the  recent  breakthrough  in  SOFC  tech¬ 
nology  and  the  investigation  of  reversible  SOFC  [4],  Most  SOEC 
studies  are  carried  out  with  the  same  cell  as  SOFC  in  order  to  benefit 
of  a  reversible  system  able  to  convert  both  hydrogen  to  electricity 
(power  generation)  and  electricity  to  hydrogen  (energy  storage) 
[5,6].  However  performance  still  remains  limited  in  electrolysis 
mode  compared  to  fuel  cell  one  [7], 

To  provide  a  better  comprehension  of  electrochemical  process, 
literature  reports  at  least  two  testing  methods.  The  first  one  con¬ 
sists  in  testing  SOEC  under  different  operating  conditions  at  a 
single  electrode.  However  it  is  impossible  to  obtain  any  polariza¬ 
tion  curve  of  the  electrode  and  to  perform  any  electrochemical 
kinetic  study.  On  the  contrary  the  second  method,  called  three- 
electrode  technique  which  involves  a  reference  electrode  located 
at  a  suitable  place,  allows  such  study  [8],  Polarization  curves  are  not 
sufficiently  relevant  to  determine  either  the  limiting  processes  or 
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Nomenclature 

avads  specific  adsorption  area  (m-1 ) 

aviPB  specific  electrochemical  area  (m-1 ) 

Cd|  double  layer  capacity  ( F  m-2 ) 

Q  concentration  of  specie  i  (mol  m-3 ) 

Qo  concentration  of  specie  i  in  the  gas  channel 
(mol  m-3) 

dg  mean  grain  diameter  (m) 

effective  molecular  diffusion  coefficient  (m2  s-1 ) 

Dk  j  Knudsen  diffusion  coefficient  (m2  s-1 ) 

dp  mean  pore  diameter  (m) 

F  Faraday  constant  (=96,485  C  mol-1) 

Jceii  three-electrode  cell  current  density  (A  cm-2 ) 
kads  adsorption  kinetic  constant  (Pa-1  s-1 ) 

kdes  desorption  kinetic  constant  (s-1 ) 

ki/_i  kinetic  constant  (m2  s-1  mol-1 ) 
kj  kinetic  constant  (s-1 ) 

L  electrode  thickness  (m) 

Mj  molecular  weight  (g  mol-1 ) 

N,  molar  flux  ( mol  s-1  m-2 ) 

rij  number  of  electrons  exchanged  at  the  jth  step 

P  pressure  (Pa) 

P,  partial  pressure  (Pa) 

R  gasconstant(=8.314JK-1  mol- 1 ) 

v;  surface  reaction  kinetic  (mol  s-1  m-2 ) 

T  temperature  (K) 

yi  molar  fraction 

a  symmetrical  factor 

e  porosity 

0cM  potential  of  the  electron  conducting  phase  (V) 

0cs  potential  of  the  ionic  conducting  phase  (V) 
rj  electrode  overpotential  (V) 

tpi  volume  fraction 

Vi  rate  of  the  zth  step  (mol  m-3 ) 

0  coverage  of  adsorbed  OH- 

02  coverage  of  adsorbed  H20 

aefi  effective  conductivity  (S  m-1 ) 

a  ionic  or  electronic  conductivity  (S  m-1 ) 

r  gas  phase  tortuosity 

rs  solid  phase  tortuosity 

r  maximal  adsorbed  specie  concentration  (mol  m-2 ) 


the  reaction  mechanism.  Thus,  another  technique,  electrochemical 
impedance  spectroscopy  (EIS)  is  used  to  identify  and  characterize 
the  various  phenomena  involved  in  electrode  reactions  in  order 
to  optimize  the  entire  process  [9,10],  Concerning  the  hydrogen 
electrode,  Ni-cermet  was  shown  to  offer  poorer  performance  in 
SOEC  mode  than  in  SOFC  mode.  Under  high  steam  concentration, 
Ni-cermet  degrades  for  low  hydrogen  content  in  the  gas  [4,7], 
However,  some  previous  investigations  performed  at  T=1000°C 
with  a  Ph2o/Ph2  ratio  =  10  demonstrated  that  nickel  electrode 
exhibits  attractive  performances  in  high  temperature  electrolysis 
operation.  Schouler  et  al.  [11]  suggest  several  possible  reaction 
schemes  for  water  electroreduction  depending  on  the  overpoten¬ 
tial  region.  For  low  overpotentials,  a  direct  electrochemical  water 
reduction  is  assumed  whereas  for  high  overpotentials  electrolyte 
redox  couple  seems  to  participate  to  water  reduction.  According 
to  Eguchi  et  al.  [4]  the  mechanism  of  steam  reduction  does  not 
only  depend  on  electrode  material;  temperature  and  moisture  par¬ 
tial  pressure  appear  to  be  critical  parameters  of  electrochemical 
reaction.  Marina  et  al.  [7]  present  Tafel’s  parameters  for  different 
experimental  conditions.  Ni-cermet  cathodic  polarizations  exhibit 


a  transport-controlled  behavior.  However  the  hydrogen  electrode 
is  less  studied  than  the  oxygen  one  because  it  is  assumed  that  the 
major  limitation  comes  from  the  anode  side  in  electrolysis  mode 
[12]. 

Computing  simulation  appears  to  be  one  of  the  most  efficient 
approaches  to  analyze  the  coupled  mechanisms  of  SOEC  opera¬ 
tion.  Modeling  of  SOFC  electrodes  could  be  focused  on  optimization 
or  characterization  [13],  It  opens  the  possibility  of  testing  the 
influence  of  parameters  such  as:  intrinsic  conductivities,  particle 
size,  graded  or  homogeneously  distributed  porosity  or  composi¬ 
tion  [14],  The  model  employed  should  depend  on  which  type  of 
loss  dominates,  on  operating  conditions  and  other  parameters  such 
as  microstructure  or  intrinsic  conductivities  [15,16],  Nevertheless, 
few  works  are  available  in  literature  regarding  the  water  reduc¬ 
tion  kinetic  on  Nickel/YSZ  cermet.  Two  electrochemical  models 
were  developed  for  an  electrolyte  supported  cell  and  simulations 
were  compared  to  experimental  data  [17,18],  Both  models  seem 
to  provide  quite  accurate  predictions  although  there  are  different. 
Only  one  theoretical  study  was  proposed  for  a  cathode-supported 
SOEC,  but  without  any  comparison  with  experimental  data  [19], 
The  major  part  of  available  models  uses  a  Butler-Volmer’s  law  for 
activation  overpotential  calculation.  However,  this  approach  seems 
to  be  not  precise  enough.  In  our  previous  work  [20],  it  has  been 
emphasized  that  a  more  accurate  electrochemical  model  for  water 
reduction  is  required  to  describe  the  high  current  density  operation 
at  low  inlet  water  concentration. 

In  this  study,  Ni-cermet  is  synthesized  and  tested  as  solid  oxide 
electrolyzer  cathode  using  three-electrode  technique.  YSZ  is  used 
as  the  electrolyte  and  Pt  as  the  counter  electrode.  The  experimental 
part  is  focused  on  the  behavior  of  cathodic  Ni-cermet  polariza¬ 
tion  in  3%H20/H2  and  3%H20-3%H2/N2  between  700  °C  and  900  “C 
with  EIS  measurement.  Polarization  curves  and  impedance  spectra 
have  been  analyzed  for  different  gas  compositions.  Indeed,  sev¬ 
eral  capacitive  contributions  have  been  recorded  depending  on 
gas  composition,  temperature  and  current  density.  In  the  model¬ 
ing  part,  a  continuum  model  is  proposed  for  a  porous  composite 
cathode  (Ni-YSZ  cermet).  It  provides  a  complete  description  of  the 
electrode  structure  and  the  processes  occurring  therein  accounting 
for  activation,  mass  and  charge  transports.  The  Ni-cermet  manu¬ 
factured  in  this  work  has  been  characterized  to  provide  the  classical 
parameters  of  microstructure  (porosity,  mean  grain  diameters) 
used  in  the  computational  approach.  Kinetic  parameter  value  of 
electrochemical  and  chemical  reactions  are  fitted  using  numerical 
optimization  method.  Predictions  from  simulations  are  compared 
to  experimental  results.  The  model  could  exhibit  the  polarization 
curves  with  one  or  two  electrochemical  or  chemical  steps  occurred 
in  water  reduction.  From  this  observation,  one  of  the  proposed 
water  reduction  processes  is  selected  to  predict  the  cell  electrical 
behavior  in  steady  state. 


2.  Experimental  part 

2.1.  Sample  preparation  and  characterization 

A  commercial  Ni  carbonate  powder  (Acros  Organics)  was  used 
as  NiO  precursor.  CH4Ni307-xH20  was  calcined  at  1000 °C  for  1  h 
to  obtain  pure  NiO.  Hydrogen  electrode  ink  was  prepared  using 
40wt.%  NiO,  60wt.%  YSZ  powders  and  the  right  amount  of  ethy¬ 
lene  glycol.  A  20  mm  diameter  and  2  mm  thickness  dense  YSZ  pellet 
(Tosoh,  8  mol.%  yttria)  was  used  as  the  electrolyte.  The  ink  was 
sprayed  over  the  substrate  using  a  3D  robot  spray  coater,  insuring  a 
good  control  of  electrode  thickness.  A  Pt  ink  was  used  to  make  the 
counter  electrode  (CE)  and  a  Pt  wire  as  the  reference  electrode  (RE). 
This  electrode  was  located  at  a  distance  close  to  half  the  electrolyte 
thickness. 
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Fig.  1.  Scanning  electron  micrograph  of  the  cross-section  of  electrolyte/Ni-YSZ  ( 


Current  density  /  mA  cnr2 


Fig.  3.  Polarization  curves  of  the  working  electrode,  Ni-YSZ  for  Ph2o/Ph2  =  1  at 
700 °C  (♦),  750 °C  (■),  800 °C  (a),  850°C  (•)  and  900°C  (+)  and  0.03  at  700°C  (0), 
750  °C  (□),  800  °C  (A),  850  ”C  (o)  and  900  °C  (-). 


Ni-Cermet  (WE) 


YSZ  Electrolyte 

1  I _ Pt  (CE) 

Fig.  2.  Schematic  view  of  the  experimental  three-electrode  cell. 


A  thickness  of  32  |jtm  for  the  Ni-YSZ  cermet  electrode  is 
estimated  from  scanning  electron  micrographs.  Image  shows 
homogenous  grain  size  throughout  the  electrode  (Fig.  1 ).  One  may 
note  presence  of  aggregates  which  confirm  that  the  electrode  is  not 
optimized  for  high  performance.  Indeed,  the  goal  of  this  study  is  to 
have  a  better  understanding  of  charge  transfer  processes. 


ratio  of  0.03.  These  results  are  in  agreement  with  Eguchi  et  al.  [4] 
studies. 

The  main  experimental  results  are  gathered  in  Table  1 .  It  appears 
that  for  low  water  content  the  electrochemical  behavior  is  clearly 
dissimilar  between  800  °C  and  900 °C.  From  700  °C  to  800  °C  the 
impedance  diagrams  involve  only  two  capacitive  contributions 
while  three  capacitive  contributions  should  be  presented  at  900  °C. 
In  addition  the  apex  frequencies  of  the  first  capacitive  contribution 
increase  from  2  kHz  to  5  kHz  with  temperature  except  at  900  °C 
where  the  apex  frequency  is  equal  to  2  kHz.  Reaction  kinetic  is 
proportional  to  the  apex  frequencies,  then  until  800  °C  the  charge 
transfer  is  enhanced  by  temperature  effect.  At  900  °C,  the  charge 
transfer  seems  to  be  inhibited  by  the  apparition  of  a  lower  charge 
transfer  process  since  there  is  more  than  one  contribution  in  low 
frequency  range.  In  addition,  observation  is  completely  different  at 
900  °C  as  shown  by  the  changing  shapes  of  impedance  spectra  at 
this  temperature  as  shown  in  Fig.  4.  The  impedance  diagram  shapes 
and  the  apex  frequencies  depend  on  steam/hydrogen  ratio  and  load. 
The  electrode  is  less  impedent  for  a  steam/hydrogen  ratio  of  1  than 
it  is  for  0.03.  The  water  effect  on  the  electrocatalytic  activity  is 
thus  clearly  pointed  out  by  this  figure.  Our  results  are  in  agree¬ 
ments  with  Schouler  et  al.  investigations  stating  that  the  transition 


2.2.  Electrochemical  characterizations 

The  setup  for  electrochemical  measurement  is  schematized 
in  Fig.  2.  The  Ni-YSZ/YSZ/Pt  cell  was  fed  with  either  3%H20/H2 
at  2.6  x  10_3Nm3h_1  or  3%H20-3%H2/N2  at  1.3  x  10“3  Nm3  h”1. 
Steam  generator  is  a  gas  bubbler  at  room  temperature.  Before  mea¬ 
suring,  the  working  electrode  is  reduced  with  dry  hydrogen  during 
24  h  at  800  °C.  All  trials  were  performed  at  atmospheric  pressure. 
The  electrodes  were  connected  to  a  Solartron  1260  for  impedance 
and  /- V measurements.  Impedance  measurements  were  performed 
scanning  from  0.1  Hz  to  65  kHz  with  applied  amplitude  of  10  mV. 

In  the  present  work,  the  three-electrode  cell  enabled  us  to  deter¬ 
mine  the  working  overpotential  at  several  points  of  the  I-V  curve. 
Impedance  measurements  were  also  carried  out  at  various  inten¬ 
sities.  These  measurements  allow  us  correcting  the  errors  arising 
from  a  shift  in  the  electrolyte  current  distribution  from  the  primary 
to  the  secondary  current  distribution  under  load  [8], 

3.  Results  and  discussion 

Fig.  3  shows  polarization  losses  in  the  Ni-YSZ  electrode.  Below 
800  °C  the  lowest  overpotentials  have  been  recorded  for  a  same 
magnitude  of  current  density  within  PH20/PH2  =  1.  At  800  °C  and 
higher  temperatures,  it  has  been  recorded  for  a  steam/hydrogen 


Table  1 

Apex  frequencies. 


T 

Ph2o/Ph2=0.03 

Ph2o/Ph2  =  1 

Jceii  =  8.6  mA  cm-2 

for  rj  =  520  mV 

Jceii  =  8.4  mA  cm-2 

for  rj  =  390  mV 

700  “C 

1st  contribution 
2nd  contribution 

2.3  kHz 

85  Hz 

1st  contribution 
2nd  contribution 

1.6  kHz 

260  Hz 

T 

Ph2o/Ph2  =  0.03 

Ph2o/Ph2  =  1 

Jceu  =  11.3mAcm- 

2  for?) =420  mV 

'  ./cell =  1 5.4  mA  cm- 

2  forq  =  360  mV 

750°C 

1st  contribution 
2nd  contribution 

3.4  kHz 

65  Hz 

1st  contribution 
2nd  contribution 

3.4  kHz 

360  Hz 

T 

Ph2o/Ph2  =  0.03 

P„2o/Ph2  =  1 

Jceii  =  16.7  mA  cm- 

2  for  q  =  297  mV  ;cell  =  16  mA  cm-2 

for  q  =  380  mV 

800  °C 

1st  contribution 
2nd  contribution 

160  Hz 

1st  contribution 

5  kHz 

600  Hz 

T 

Ph2o/Ph2=0.03 

Ph2o/Ph2  =  1 

Jceii  =  75.4  mA  cm- 

2  for  77  =  300  mV  Jce\\  =  64.7  mA  cm- 

2  for  rj  =  600  mV 

900  “C 

1st  contribution 
2nd  contribution 

2  kHz 

650  Hz 

7  Hz 

1st  contribution 
2nd  contribution 

2  kHz 

500  Hz 

3.5  Hz 
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Fig.  4.  Impedance  spectra  for  a  working  electrode  potential  of  -100  mV  and  -1  V/ref  with  steam/hydrogen  ratio  of  0.03  and  1  at  900  °C. 


between  the  varying  mechanisms  depends  on  operating  temper¬ 
ature,  gas  composition  and  applied  potential  [11],  At  900  °C,  for  a 
difference  of  - 1  V  between  the  working  and  the  reference  electrode 
(-1  V/ref),  the  low  frequency  contributions  have  the  same  ampli¬ 
tude  as  high  frequency  contributions  for  both  ratios.  First  the  low 
frequency  contribution  could  be  assign  to  a  diffusion  process.  But 
this  assumption  is  not  in  agreement  with  the  small  current  densi¬ 
ties  recorded  in  these  experiments  and  the  absence  of  nitrogen  for 
Ph2o /Ph2  =  0.03  does  not  involves  a  lower  amplitude  of  this  capaci¬ 
tive  contributions.  A  simple  diffusion  process  could  not  explain  this 
behavior.  These  results  emphasize  the  possible  existence  of  coupled 
processes. 

3.1.  Theoretical  assumption  of  reduction  process  from  literature 

According  to  literature,  two  kinds  of  reaction  mechanisms  have 
been  suggested  on  Ni-YSZ  interfaces  [1 1,21  ].  For  a  nickel  electrode 
at  low  Ph2o  the  mechanism  could  correspond  to  the  following  reac¬ 
tions: 

H20gas  +  02-B^20H-ds  (1) 

Vo  +OHa-ds  +  e-^Og  +  (1/2)H2  (2) 

and  the  corresponding  kinetics  equation  could  be  expressed  as  fol¬ 


lowing  [22]: 

V!  =  ki  r2exp  (^(0cM  -  0cs))  e2  -  fe_,r2 

(3) 

y2  =  k2rexp(^(0cM-0cs))  7  (1  -  6>)  -  Jc_2  T 

x  exp  (-(1  ~“)F(0cM  -  0cs ))  0  (4) 

At  high  steam  pressure,  a  direct  reduction  of  steam  has  been 
demonstrated.  The  mechanism  could  then  be  expressed  as: 

H2Ogas  -*•  H2Oads  (5) 

Vq  +H20ads  +  2e-^C>S  +  H2  (6) 

The  kinetics  of  steam  adsorption  and  reduction  being  similar  to 
oxygen  kinetics  in  SOFC  mode  [13]  may  be  described  by: 

V3  =  r 2kadsPH2o(l  -  02 )  -  F ikdes02  (7) 


V4  =  k4r2  exp  (“J(0cM  -  0cs))  -  e2)  -  k_4r2 

X  exp  (-(1  ~fF(0cM  -  0cs))  02  (8) 

The  rates  of  steps  (v,)  depend  on  hydrogen  or  steam  concen¬ 
tration  and  on  coverage  of  adsorbed  species  (0  for  OH-  or  02  for 
H20)  with  kinetic  constants  ( k_y ).  According  to  Schouler  et  al.,  the 
transition  between  the  two  mechanisms  appears  close  to  -  1 .2  V/ air 
for  Ni-YSZ  [11],  Moreover,  this  transition  also  depends  on  operat¬ 
ing  conditions  (temperature,  gas  composition)  such  as  observed  in 
Table  1.  The  presented  experimental  results  do  not  permit  to  select 
a  mechanism  as  a  function  of  temperature  or  Ph2o/Ph2  ratio.  Klein 
et  al.  have  emphasized  that  the  overpotential  are  not  homogenous 
within  the  porous  electrode  due  to  the  material  electric  conductiv¬ 
ity  [23],  Then,  it  is  quite  possible  that  both  electrochemical  ways 

In  other  hand,  some  limitations  due  to  an  electronic  conductivity 
of  the  YSZ  should  appear  in  SOEC  mode.  According  to  Pham  and 
Glass  [24],  an  electronic  conduction  due  to  partial  electrolysis  of 
the  electrolyte  should  occur  for  an  oxygen  pump.  However,  this 
effect  at  “high  enough  polarization”  (e.g.  10-1 5  atm  at  800  °C)  is  still 
much  higher  than  values  reported  by  Yuan  and  Kroger  (10-27  atm) 
[25].  Subsequently,  in  the  context  of  SOEC,  it  was  reported  at  high 
concentration  polarization  losses,  possible  electronic  short-circuit 
within  the  electrolyte  [26], 

Although  in  SOEC  mode,  the  cermet  assumes  a  cathodic  polar¬ 
ization  which  protects  the  nickel  material.  Accorsi  and  Bergmann 
[27]  have  suggested  that  steam  partial  pressure  in  SOEC  mode  must 
be  less  than  20%  to  avoid  the  nickel  oxidation.  And,  it  is  established 
that  the  nickel  could  be  not  oxidized  if  the  H2  composition  is  higher 
than  2%  at  900  °C  [28],  The  possible  re-oxidation  process  of  cermet 
involves  a  degradation  of  the  hydrogen  electrode  [  29  ]  and  the  risk  of 
Ni  re-oxidation  increases  with  decreasing  Ni  grain  size,  the  gradient 
of  overpotential  and  steam  composition  through  the  electrode.  So, 
the  possible  reduction  of  YSZ  and  re-oxidation  of  Ni  should  impact 
the  electrochemical  process  of  steam  reduction. 

4.  Modeling  part 

4.1.  Continuum  model 

The  continuum  model  [13]  is  isothermal  and  one-dimensional 
(Fig.  5).  It  is  assumed  that  the  gas  partial  pressure,  concentration 
of  gaseous  and  adsorbed  species  and  overpotential  are  uniform 
over  the  geometrical  surface  area  of  the  electrode.  Ohmic  losses  in 
the  well-connected  electronic  conductor  phase  are  neglected.  Thus, 
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this  phase  is  regarded  as  equipotential.  Transport  of  oxygen  vacan¬ 
cies  in  YSZ  is  strictly  ohmic.  In  the  case  under  study,  water  steam 
is  provided  with  hydrogen  and  in  some  conditions  with  nitrogen. 
Stefan-Maxwell  diffusion  model  [30]  is  used  regarding  the  signif¬ 
icant  difference  between  molecular  weight  of  species.  Moreover, 
Knudsen  diffusion  [31]  should  be  considered  in  the  electrode  that 
is  a  porous  media.  Thus,  the  Dusty-Gas  Model  (DGM)  is  used  at  cath¬ 
ode.  Such  model  seems  to  provide  the  most  accurate  description  of 
mass  transport  for  SOFC  anode  [32], 

Molar  fraction  variation  of  specie  i  y,(x,  t)  is  calculated  solving 
the  following  equation: 


to 

dx 


O) 


The  Knudsen  diffusion  coefficient  of  the  species  z  (Dfc  i)  is  cal¬ 
culated  from  relation  (Eq.  (10))  where  the  mean  pore  diameter  is 
determined  considering  electrode  characteristics  (Eq.  (11)). 


Dk 

dp 


/  8  RT 

V  nMi 


2  e 
3l^£ 


dg 


(10) 

(11) 


The  effective  diffusion  coefficient  is  the  diffusion  coefficient  in 
a  porous  media.  Different  models  were  used  to  assess  this  coef¬ 
ficient  [33,34],  Our  models  developed  for  SOEC  cathode  use  the 
Bruggeman’s  law. 


Du  =  Du£T 


(12) 


The  overpotential  (rj)  is  equal  to  the  difference  between  poten¬ 
tial  of  the  electron  conducting  phase  and  the  potential  of  the  ionic 
conducting  phase.  Then,  continuity  equation  related  to  faradic  pro¬ 
cess  and  double-layer  charging  process  is  given  by: 

ai>rPBCdi  ^  -  9  ^  +  avTPBF^Tijrj  =  0  (18) 

j 

apppB  is  specific  electrochemical  area  which  represents  elec- 
trochemically  active  surface  per  unit  volume  of  electrode.  This 
area  could  be  estimated  from  statistic  law  of  binary  mixture  pow¬ 
der  [36,37],  Assuming  a  one  rate-limiting  step,  the  Butler-Volmer 
law  can  be  used  to  describe  the  kinetic  as  demonstrated  for  the 
Volmer-Heyrovsky  mechanism  [38],  A  parametric  study  revealed 
a  similar  behavior  for  the  global  and  the  approximated  kinetics. 
For  calibration  and  validation  steps,  the  electrochemical  kinetic  has 
been  described  by  the  following  equation: 

(19> 

Electrode  is  a  composite  of  electronic  and  ionic  conductors. 
Effective  electronic  and  ionic  conductivities  are  estimated  by  the 
following  relation  [39]: 

aef  or  io  =  ^Ni  or  YSZ~ or  YSZ  (20) 

Ionic  conductivity  of  the  YSZ  is  considered  as  temperature 
dependent  by  the  equation  [40]: 

ovsz  =  0.334  x  105  exp  (-10’^00)  (21) 

Nickel  electronic  conductivity  is  taken  as  a  constant  since  its 
variation  on  the  studied  temperature  range  is  negligible.  A  value  of 
6.5  x  104  S  m-1  has  been  considered  [41  ].  A  potential  is  set  at  the 
electrode/electrolyte  interface  and  Ionic  current  is  nil  at  the  gas 
channel/electrode  interface  (x  =  0). 

V(x=L)  =  rf  (22) 


Mass  conservation  for  each  species  could  be  written  as  follows: 


9Nh2  ,  v-  P  dyu2 

dx  ±aVTPBZ^vi-  RT  dt 

(13) 

3Nh2o  P  9Yh2o 

dx  ±wtPB&=RT  dt 

(14) 

9%2  P  9yN2 
dx  ~  RT  dt 

(15) 

Total  pressure  is  assumed  to  be  constant  then  one 
the  following  condition: 

can  consider 

2>-* 

(16) 

In  steady-state,  mass  balances  are  solved  using  Runge-Kutta 
method  and  charge  balance  using  Newton-Raphson  in-house 
developed. 

4.2.  Simulation  results 

The  continuum  model  has  been  solved  to  generate  polarization 
curves.  First,  kinetics  parameters  have  been  determined  in  a  cal¬ 
ibration  step.  During  this  step,  the  least  square  distance  between 
the  model  and  the  experimental  data  are  minimized.  This  para¬ 
metric  optimization  has  been  carried  out  using  the  simplex  method 


Molar  fractions  are  known  at  x  =  0.  At  x  =  L,  flux  are  nil.  Diffusion 
of  adsorbed  species,  OH_a(js  et  H20a<jS  are  assumed  to  be  negligi¬ 
ble  compared  to  diffusion  process  in  the  gas  phase  since  diffusion 
coefficients  of  the  adsorbed  species  are  on  the  order  of  1 0-12  m2  s-1 
[35].  Thus,  mass  balance  becomes: 

(17) 


Table  2 

Parameter  value  used  for  the  calibration  step. 


Parameter  Value  Unit  Reference 


TS  1.7  -  [13] 

r  4.8  -  [20] 

dg  1  x  10-6  m  This  study 

e  03  -  This  study 

tp ysz  0.6  -  This  study 

a  0.5  -  [7] 
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I  from  the  calibration  step. 


Ph2  O  /Ph2  OVrpnfc2bis  OVjpB  fc— 2bis  CIVypR^  G2  UVTPtik-4  F2  UVa^sF2l<ads  aVdds  LYfcdeS 

0.03  79.2  2.38  7.47  x  106  770  134  xlO5  4.03  x  107 

1  239  239  1.22  x  107  2.17  x  106  2.16  x  105  8.89  x  107 


contribution.  The  kinetics  of  steam  adsorption  Wads  ^2  kads  •  steam 
desorption  qvads P2  kdes  and  reduction  are  restrained  in  the  same 
magnitude  values  under  both  studied  atmospheres.  Nevertheless, 
it  is  worthy  to  note  the  significant  scale  change  of  the  fitted  value 
for  the  oxidation  kinetics  ai7rPBk_4  r2  from  770  mol  s-1  m-3  at 
Ph2o/Ph2  =  0.03  to  2.17  x  106  mols-1  m-3  at  Ph2o/Ph2  =  1.  Then, 
our  investigations  are  in  agreement  with  the  observations  of 
Sukeshini  et  al.  who  highlighted  a  strong  difference  of  electro¬ 
chemical  behavior  between  anodic  and  cathodic  polarizations  of 
Ni-YSZ  cermet  [42],  The  definition  of  the  transfer  current  density 
in  macro-scale  model  should  be  reconsidered  in  the  light  of  this 
result.  Another  equation  should  be  defined  taking  into  account  for 
example  an  adsorption  step. 

5.  Conclusion 


Fig.  6.  Comparison  between  experimental  (a)  and  simulated  polarization  curves 
considering  an  adsorbed  water  molecule  ( — )  or  an  adsorbed  hydroxide  ion  (— )  in 
the  reaction  path  under  a  Ph2o/Ph2  of  0.03  at  800  °C. 


available  in  Matlab®.  In  order  to  avoid  a  too  large  number  of  param¬ 
eter,  some  of  them  are  estimated.  Considered  values  are  gathered 
in  Table  2. 

Kinetics  parameters  determined  with  the  simplex  method  are 
given  in  Table  3.  The  best  fit  is  obtained  for  the  model  where 
an  adsorbed  water  molecule  is  involved  in  the  reaction  path  (see 
Figs.  6  and  7).  The  other  reaction  path  (Eqs.  (1)  and  (2))  has  been 
modeled  using  the  Butler-Volmer  equation  assuming  one  rate- 
limiting  step.  The  disagreement  observed  in  Figs.  6  and  7  leads 
to  the  conclusion  that  this  equation  is  not  adapted  to  model  the 
electrochemical  behavior  of  the  electrode  for  the  considered  gas 
compositions. 

Comparison  between  simulation  and  experimental  data  leads 
to  the  conclusion  that  there  is  more  than  one  rate-limiting  step 
reaction  for  water  reduction.  The  cathodic  reaction  takes  place 
at  the  triple  phase  boundaries  described  by  the  m/rpB  coefficient 
of  specific  electrochemical  area.  In  this  study  the  in-house  cer¬ 
met  exhibits  large  aggregates  which  involve  small  specific  reaction 
areas  av tpb  and  Wads  ^2  ^des  •  These  small  performances  involve  a 
high  magnitude  of  ASR  and  make  easy  the  identifications  of  varying 


In  this  study,  a  Ni-cermet  electrode  has  been  electrochemically 
characterized  using  the  three-electrode  technique.  Polarization 
curves  and  impedances  measurements  have  been  carried  out 
for  temperature  between  700  °C  and  900  °C  for  two  gas  com¬ 
positions.  The  influence  of  water  activity  on  electrode  behavior 
has  been  pointed  out.  The  electrochemical  behavior  of  the 
Ni-cermet  depends  on  temperature  and  steam/hydrogen  ratio. 
For  temperatures  below  900  °C  and  steam  content  close  to 
3vol.%,  the  polarization  curves  exhibit  lower  overpotential  with 
a  steam/hydrogen  of  1  than  0.03.  Moreover  the  electrochemical 
impedance  spectra  show  a  low  frequency  contribution,  but  the  dif¬ 
fusion  processes  which  are  classically  attributed  at  this  range  of 
frequency  are  not  valid  due  to  low  current  densities  recorded.  In 
the  second  part  of  this  work,  a  continuum  model  has  been  devel¬ 
oped.  Two  different  mechanisms  have  been  studied  for  the  mass 
and  charge  transfer.  The  first  one  has  been  modeled  assuming  one 
rate-limiting  step.  This  assumption  used  to  justify  Butler-Volmer 
equation  does  not  permit  relevant  computing  predictions  of  polar¬ 
ization  curves.  Therefore,  the  Butler-Volmer  description  is  not  valid 
to  model  water  reduction  on  Ni-cermet  for  this  water  content.  The 
mechanism  considering  the  adsorption  of  water  molecule  seems  to 
be  the  most  plausible.  When  two  rate-limiting  steps  are  modeled, 
comparisons  between  simulated  polarization  curves  and  experi¬ 
mental  ones  have  a  good  agreement. 
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